The stochastic process of equilibrium fluctuations, of a system with long range interactions 
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The relaxation towards equilibrium of systems with long range interactions is not yet understood. 
As a step towards such a comprehension, we propose the study of the dynamical equilibrium fluctua- 
tions in a model system with long range interaction. We compute analytically, from the microscopic 
dynamics, the autocorrelation function of the order parameter. From this result, we derive analyt- 
ically a Fokker Planck equation which describes the stochastic process of the impulsion of a single 
particle in an equilibrium bath. The diffusion coefficient is explicitly computed. 

PACS numbers: 



A number of physical systems are governed by long 
range interactions. Some examples are given by self- 
gravitating systems, two dimensional incompressible, or 
geophysical flows, some models in plasma physics. For 
such Hamiltonian systems, the non additivity of the in- 
teractions makes the usual thermodynamic limit N — > 
co,V ^ co irrelevant. Microcanonical average is how- 
ever still relevant, and generically leads to a mean field 
description of the equilibrium, exact in the limit N — > oo 
PJ . The relaxation toward equilibrium of these systems 
still has to be completely understood. The phenomenol- 
ogy of the dynamics shows that a rapid relaxation leads 
to the formation of quasi-stationary structures, which 
may be out of equilibrium statesfsee Q for astrophysical 
and geophysical examples and [j] for spin system ones). 
In most of cases, this is explained by the existence of sta- 
ble stationary states of the associated Vlasov equation, 
which describes the dynamics by approximating the po- 
tential by a mean field one. In such stable situations, the 
Vlasov dynamics is a good approximation of the particle 
dynamics, on typical time scales diverging with the num- 
ber of particle |j| . The relaxation towards equilibrium of 
these structures is then associated to the fluctuations of 
the potential around its mean field approximation, and is 
thus very slow. One of our goal is to understand such a 
relaxation, which is of particular interest, for instance in 
the study of astrophysical structures, turbulence parame- 
terization in geophysical flows, etc. Some works towards 
a kinetic description of this relaxation have been pro- 
posed, for instance, by Chandrasekar in the context of 
self-gravitating systems j3|, Chavanis for the point vor- 
tex model Q or for the two-dimensional Euler equation 
0, or in plasma physics Q. In each of these cases, the 
relaxation is then described by a Fokker-Planck equa- 
tion or some generalizations. These results are mainly 
obtained by formal considerations and the diffusion coef- 
ficient in the Fokker-Planck equation is always expressed 
as a Kubo formulae or, equivalently, by a formal integral 
of the Liouvillian of the dynamics. The diffusion coef- 
ficient has been computed, in some limits, for the point 
vortex model Q and for self- gravitating systems. 

In the kinetic theory of dilute gases, the Boltzmann 



equation has lead to the computation of transport coeffi- 
cients Q . This is an example of explicit computation of a 
diffusion coefficient for a system, with a large number of 
particles. A complete mathematical proof of this result 
directly from the Hamiltonian dynamics is however still 
to be achieved. The computation of the diffusion coeffi- 
cient for the standard map is a classical example, for 
a system with a small number of degrees of freedom. On 
the past decades, the issue of the link between chaotic 
Hamiltonian dynamics and diffusive properties has been 
addressed on a general framework ]. We also note 
works on the relaxation to equilibrium of a massive piston 
in interaction with two out of equilibrium perfect gases 
[ri| . which a Vlasov like behavior. 

We will show that the diffusion coefficient for systems 
with long range interactions can be computed in the large 
density limit (N — > oo with a fixed volume and renormal- 
ized interaction). At statistical equilibrium, one obtains 
the mean field description typical for long range inter- 
acting systems. Near the equilibrium, particles have an 
integrable motion, perturbed by the fluctuations of the 
mean field around its equilibrium value. This leads to 
the relaxation towards equilibrium. The self consistent 
nature of the fluctuations (the mean field oscillates due 
to small particle deviations, themselves due to the mean 
field fluctuations) is however an essential feature of this 
process. 

In order to precise these ideas, we consider a simple 
toy model of long-range interacting system: the Hamil- 
tonian Mean Field model (HMF). In this framework, as 
a first step towards the study of the relaxation towards 
equilibrium, we consider the equilibrium dynamical 
fluctuations. We first propose an analytic computa- 
tion of the autocorrelation function of the mean field 
order parameter. From this result, we can derive a 
Fokker-Planck equation which describes the stochastic 
process of a particle in interaction with a bath of N — 1 
particles in equilibrium. The diffusion coefficient is then 
explicitly computed, from the microscopic dynamics. 
We finally conclude by discussing generalization to out 
of equilibrium situations, and more realistic models. 
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The Hamiltonian of the attractive HMF model 121 is 
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Because of its simplicity, a large number of authors 
have considered this model and its repulsive counter- 
part (with the opposite sign for the potential energy). 
The HMF model is the "harmonic oscillator" for long 
range interacting systems. We refer to [Isj for a re- 
view. Let us define the magnetization M by 7VM = 
Y,k=i eie " ( M = M * + iM v)- Because the kinetic en- 
ergy per particle e c may be exactly expressed as 2e c = 
2E — 1 + M 2 (E is the energy per particle), and be- 
cause M is a simple sum of N variables, the computa- 
tion of the static microcanonical quantities is straight- 
forward. For instance we obtain the volume of the 
shell of the phase space, with energy E: ft(E) ocjv— >-oo 
/p 1 dM B (M) exp (NS (E, M)) with the entropy S given 
by S (E, M) = C (M) + log (2E + M 2 - l) /2, where 
C(M) = log(I (V>(M))) - Mip(M), I is defined by 
2nIa(M) = Jq* d6 exp (M cos 0), and tp as the inverse 
function of dloglo/dM. The use of the saddle point 
method, in the previous integral, shows that an over- 
whelming number of configurations have a magnetiza- 
tion close to the equilibrium value M e (E) defined by 
dS (E, M e ) /dM = 0. This equation shows that, above 
the critical energy E c = 3/4, whereas below E c a sec- 
ond order phase transition occurs. The density in the 
fi— phase space (all angles and momenta are projected 
on a {0,p) space) may be evaluated as Je (p, 0) oc 
exp (— (3 (E) (p 2 /2 — M e (E) cos#)) , up to a translation 
of angles. These results are equivalent to the canonical 
ones (see [l2T|l. 

In the following we consider only energies greater than 
the critical one E > E c . In such a case, the equilibrium 
is homogeneous : M e = 0. Then (3 = 1/ (2E - 1) and 
f(p,0) oc exp (— Pp 2 /2). The static fluctuations of M 
may also be computed, from the second derivative of S 
with respect to M, at the equilibrium point. We obtain 
a Gaussian magnetization with N (M*M) = 2/ (2 - 0) 
and (M*M*) = (MM) = (M* is the complex conju- 
gate of M). The magnetization has typical fluctuations 
of order iV -1 / 2 , we thus re-scale it accordingly, by defin- 
ing : 
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The aim of this letter is to study the dynamical equilib- 
rium fluctuations of this system. From the Hamiltonian 
(QJ, one obtains the equations of motion: 



d0 k , dp k N~ 1 / 2 

_ =pk and _ = 



im* (t) + CC) (3) 



(CC means the complex conjugate of the previous ex- 
pression). From the motion equation (pj, thanks to the 
smallness of the mean field fluctuations, the motion of 
any particle may be treated perturbatively in the limit 
N — > 0. We expand the variables in power of iV -1 / 2 : 
k = e kfl + N~ l / 2 6 kA + Pk = Pk,o + N'^-^pk.i + ... 
and m = m + N~ 1 / 2 vcii + ... (the magnetization l|2| has 
to be self-consistent) . The zero order motion is a free bal- 
listic one : p k ,o (t) = pi and 6k,o (t) = 0% +p a k t, where 9% 
and p° are the values of 6 and p for t = 0. The expression 
<EJ clearly shows, that such a perturbative description, 
around this simple zero order dynamics, will remain valid 
as soon as t <§; N 1 / 2 . This expansion leads, to the first 
order, to Ok.i (t) = J dupk t i (u) with : 



PkA (t) 



du 



i* (u) e« +p ° u ) + CC 



(4) 



A peculiarity of this asymptotic expansion, is that the 
magnetization m is a sum of N variables, where 

N^ 1 / 2 is the expansion parameter. A sum of N order 
N^ 1 / 2 terms may be of order 1. To obtain the zero or- 
der magnetization, we thus have to include the first order 
expression of the angles 9k,i- We then obtain: 

mo (t) = a — 1/2 I du J dv [6hiq (v) — cmo (v)] (5) 

J JO 

where a(t,6° k ,p° k ) = N~^ 2 eW+X*) , 

b(t,v,6° k ,pl) = N -i E ^ ie KK+ P l(t +v)) and 

c(t,v,6l,p%) = iV- 1 X)f =1 e i, * ( * _, ' ) . The expression 
® clearly reflects the self-consistent nature of the 
motion : the magnetization at time t depends on the 
magnetization at previous times. 

From Ipjfl. we compute the magnetization autocorre- 
lation function at leading order (hiq (t) mo (0)), where 
the bracket denotes microcanonical averages on the vari- 
ables (8 k ,p ( l)- We first note that b and c are equal to 
their microcanonical average plus fluctuations of order 
^y-i/2 These fluctuations can be neglected at the or- 
der considered, b and c are thus treated as independent 
of the magnetization. Some lengthy computations lead 
to (a(t,e°, P °)mo) =2cxp(-t 2 /2/3)/(2-/3), (b(t,v, 0°,p°)) = 
and (c(t,v,e° k ,p° k )) = 2e x p(-(t-t>) 2 /2/3) /(2-/3) (up 

to order N^ 1 / 2 corrections). From (jjj, we then ob- 
tain (m*(i)m(0)) = 2<t> (t) / (2-/3) and <m*(t)m*(0)> = 
(m (t) m(0)) = 0; where the function <p is given by the so- 
lution of the integral equation : 



. (t) = exp - 



2/3 



dv v exp 



2/3 



4>{t-v) 



(6) 

We remark that the rhs integral is a convolution. This 
makes the solution of this equation by a Laplace trans- 
form natural. We do not report the result. Whereas the 
first term on the rhs of this integral equation is due to 
the integrable ballistic motion of the particles, the second 
term reflects the self-consistent nature of the dynamics. 



3 



To our knowledge, it is the first derivation of such a mem- 
ory term, directly from the dynamics, in an Hamiltonian 
system with a large number of particles. 

To have a physical insight on this autocorrelation func- 
tion, we compute the asymptotic behavior of <j>. Firstly 
4> (t) oq-,0 ex P ( _ (2 — p) t 2 / (4/3)). This approximation 
is obtained using the value of the second derivative of <j> 
in : 4>i. FromUjJ by a Taylor expansion, we compute 
4>x = (2 — 0) I (2/3). Such a Gaussian behavior for small 
times, would be typical of a ballistic behavior. However, 
we note that the coefficient (2 — /3) / (2/3) is not uniquely 
due to the integrable zero order dynamics, but is renor- 
malized by the memory term. Secondly we obtain : 

(t) oc*^ A 09) exp (-7 (/?) t) ; 7 (0) = (2//3) 1/2 F" 1 (/?) , 

(7) 

where F 1 is the inverse of the function F, with F (x) = 
2/ (l + y/irx exp (x 2 ) erfc (— x)) , where erfc is the com- 
plementary error function This exponential limit for the 
autocorrelation function is natural, as it corresponds to 
the Markovian limit for the magnetization stochastic pro- 
cess. Using the expression @ as an Ansatz and eval- 
uating the integral equation © at dominant order for 
t — ► oo, it is possible to obtain this result for 7. The 
lower inset of Fig. shows the relaxation constant 7 
as a function of (3. Near the critical energy (/3 C = 2), 
the relaxation constant tends to 0. This indicates that 
near the critical point, the relaxation time diverges. On 
the contrary, for large energy, 7 diverges and the relax- 
ation time is very small. Let us compare these results 
for the autocorrelation function with numerical results. 
We first numerically compute the theoretical prediction 
for 4> (0 (t) may be numerically computed either from its 
Fourier transform, its Laplace transform or by an itera- 
tive scheme directly from Fig. Qshows a comparison 
of this theoretical autocorrelation function with the one 
obtained directly from the integration of the Hamiltonian 
dynamics JQ. 

Because the stochastic process is stationary, using the 
Wiener-Kinchin theorem, the spectral density of the com- 
plex magnetization may be computed from the Fourier 
transform of the autocorrelation function. As the integral 
equation describing the autocorrelation function is a con- 
volution, the computation of this spectral density is easy. 
Defining S (u>) = 2/n J °° dt cos (wt) (m* (t) mmO)), one 
obtains : 

= 8(^/2) 1 / 2 exp(- / 3^ 2 /2) 

2 [(2 -/3) + P 3 / 2 loA {P 1 / 2 ^)] 2 + tt/3 3 w 2 exp {-/3lo 2 ) 

(8) 

where A (x) = exp (-x 2 /2) J* du exp (u 2 /2). We have 

S(u) ~„->oo (2/3/7r) 1/2 exp {-(3uj 2 /2). 

Let us now consider the diffusion of the momentum 
p of a single particle, where all others particles have a 
random angle and momentum according to the micro- 
canonical distribution (one particle in a bath at equilib- 
rium) . Let us denote ( Ap) (p, t) the mean displacement 
p(t) — p(0), and (Ap 2 ) (p, t) the mean square displace- 
ment of a particle knowing that its initial momentum 
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Figure 1: The magnetization autocorrelation function: the 
predicted value 10 and the numerically computed value are 
both represented. They are indistinguishable (maximum ab- 
solute error of 3.10" 3 ). We have used E = 2.5, f3 = 
1/ (2E — 1), N = 10 000, and averaged over 18 samples, each 
one of duration t — 8000. The upper inset shows that the 
exponential decay of the autocorrelation function is a good 
approximation times greater than 2 or 3, for these param- 
eters. The lower inset shows the relaxation constant as a 
function of the inverse temperature 7 {0) (sec0. 

(p (0) = p) . From |@J and the results for the autocorre- 
lation function or for similar quantities, it is possible to 
compute explicitly (Ap 2 ) (p, t), at the leading order in N, 
for any time such that t <C N 1 ^ 2 (perturbative descrip- 
tion of the dynamics). The quantity (Ap 2 ) (p,t) has a 
transient behavior on a time scale of order 1 (the explicit 
computation is feasible, but not reported), followed for 
1 *C t <C N 1 / 2 , by a diffusive behavior. We then obtain 
(Ap 2 ) (p,i) ~ 1<<t « w -i/2 2D(p)N-H, with 

1 r 00 

D (p) = o / dt ( m * (*) m (°)) cos (P*) ( 9 ) 

2 Jo 

This result is the equivalent of a Kubo formulae. However 
it states a bite more: the diffusion coefficient is there ex- 
pressed as the autocorrelation of the mean field and not 
as the autocorrelation of the force. We note that the dif- 
fusion coefficient is proportional to the spectral density : 
D (p) = irS (p) /4. This is a peculiarity of this model for 
which the interaction is built with a cosine. An analytical 
expression for D is thus obtained from {HJ. The compu- 
tation of (Ap) (p, t) may be done following the same pro- 
cedure. Please note, however, that the TV -1 / 2 contribu- 
tion vanishes. The lower order contribution comes from 
a perturbative description of the dynamics at order 2 
(pfc,2 {t) and mi (t)). The systematic momentum change 
is N (Ap) (p, t) ~ 1<t<<JV -i/2 {dD (p) I dp - ppD (p)) t. In 
the following, we will see that this result can be deduced 
from the stationarity of the stochastic process (or equiv- 
alent^ from the fact that the distribution for p tends 
to the equilibrium distribution). This is the equivalent 
of an Einstein relation. Fig. El shows that the analyt- 
ical diffusion coefficient with the numerically computed 
N (Ap 2 ) (p, t) ft agree. 

We have observed a diffusive behavior for the momenta 
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Figure 2: The solid curve shows the mean square displacement 
of a particle in function of its initial momentum , normalized 
by N and divided by the time (N ( Ap 2 ) (p, t) / (2t)) (this is 
not a distribution), for four values of time : t — 10, 15, 20, 25 
; N = 10 000, f3 = 1/4 . As the curves are superposed for 
time, this shows that the motion is actually diffusive. The 
dashed curve represents the predicted result D (p) (equation 
ISlwith D (p) — ttS (p) /4), no fit. This confirms the theoretical 
analysis, up to errors due to an incomplete statistics. 

with a systematic momentum drift, for 1 <C t <C 
TV 1 / 2 . Moreover the mean displacement and the mean 
square displacement are small as they scale like TV -1 . 
These two facts are the two hypothesis for the derivation 
of a Fokker-Planck equation. Thus any momentum dis- 
tribution function / (p) evolves, at the leading order in 
N, through the equation: 

f=^IH(f+<w)) (io) 

This equation is valid for time t 1. For the deriva- 
tion of the mean square and mean displacement, we have 
assumed t <C N 1 ^ 2 (perturbative description). However, 
the previous analysis has also shown that the correla- 
tion function decays exponentially for large time. The 
correlation time for the force (or equivalently the mag- 
netization) is then of order 1 and is thus much smaller 
than N 1 / 2 . This is a first indication that the stochastic 
process may becomes Markovian for time much smaller 
than N 1 / 2 . If it is actually such, the Fokker-Planck 
will be correct for any time t. We note that this equa- 
tion actually converges towards the equilibrium density 
P eq (p) = (/3/2^) 1/2 exp(-/3p 2 /2). 

In this letter, we have analytically computed the auto- 
correlation function for the HMF model. We have used 
this result to compute analytically the diffusion of the 
momentum of a single particle in an equilibrium distri- 
bution. We have obtained a Fokker-Planck equation for 
which the diffusion coefficient is explicitly computed. A 
more complete study of the magnetization stochastic pro- 
cess, the detailed computations, and the study of this 
Fokker-Planck equation will be addressed in a forthcom- 
ing paper [14]. Due to the asymptotic decay of the dif- 
fusion coefficient, for large momentum, the spectrum of 
the linear operators of the Fokker-Planck equation has no 
gap between the eigenvalue corresponding to the ground 
state and the other eigenvalues. In the introduction, we 
have introduced this work as a first step towards the de- 
scription of relaxation towards the statistical equilibrium. 



We hope to generalize in the future these results to such 
out-of equilibrium states. Let us note however, that the 
Fokker-Planck equation should correctly describe the re- 
laxation towards equilibrium of states sufficiently close 
to equilibrium. The generalization of this letter results 
to other long range interacting particle models may fol- 
low the same path, with technical difficulties linked to 
the continuous nature of the mean field in most systems, 
and the theoretical problems linked with the divergence 
of some interactions at small scales (point- vortices, self- 
gravitating systems and plasma). 
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pean network Stirring and mixing, RTN2-2001-00285. 
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